PROCEEDINGS 

OF SCIENCE 



oo 

o 
o 

(N 
•*-> 

o 

o 



> 



oo 

o 

> 

• 1—1 

X 



Block renormalization group transformations and 
overlap fermions 



Nigel Cundy 

Universitdt Regensburg, Universitat Strafle 31, 93040 Regensburg 



E-mail: nigel . cundy@physik . uni-regensburg . de 



In this preliminary work, I provide the outline of an argument (leaving the full proof to a future 
publication) that there exists a valid renormalization group blocking transformation which con- 
verts the continuum fermion action into a Ginsparg-Wilson lattice action. I construct the blocking 
for the massless overlap operator as a specific example, indicating how other Ginsparg-Wilson 
lattice Dirac operators can be derived in a similar fashion. This renormalization group transfor- 
mation modifies the gauge action and adds a number of irrelevant terms to the lattice action. The 
procedure is not valid for lattice Dirac operators which do not exactly satisfy the Ginsparg-Wilson 
relation, for example the Wilson operator. 
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1. Introduction 

Since chiral symmetry determines many of the low energy properties of QCD, lattice simu- 
lations should ideally respect chiral symmetry to a high accuracy. The Ginsparg- Wilson frame- 
work [jl|, ||] offers a method of constructing lattice Dirac operators with an exact 'chiral symmetry.' 
Examples of Ginsparg-Wilson lattice Dirac operators have been known for some time, and the 
overlap operator, which exactly satisfies the Ginsparg-Wilson chiral symmetry, is currently being 
used in lattice simulations of two flavour ^ and 2+1 flavour [Q] lattice QCD. 

However, it is long been known that the overlap []|] and perfect action [§, [7|] fermions are not 
the only possible Ginsparg-Wilson Dirac operators, and there is some theoretical interest in trying 
to describe the connection between the various possible solutions to the Ginsparg-Wilson equation. 
In [0], I described the derivation of a group of valid chiral lattice Dirac operators, based on the 
overlap operator, all of which satisfy a Ginsparg-Wilson relation. The general form for massless 
Dirac operators of this type was given as 



D' N 



*(75£) 



h q 



1+r 



^(Y5£ + £Y5) 



Ys£ , (1-1) 



where s, t, h, q and r are analytic functions, only constrained by the need to preserve the correct 
continuum limit and give the doublers an infinite mass, and e = sign (ysD^), where Dy/ is a kernel 
operator, and, as with overlap fermions, can be any valid lattice Dirac operator at a negative mass 
(for example the Wilson operator). 

It remains unlikely that any of these operators offer advantages over standard overlap fermions. 
However, there is one intriguing case, where r[x] = l/\x\ (neglecting that proof of the locality of 
the Dirac operator given in [§] depends on the analyticity of r; although that the proof breaks 
down does not necessarily mean that the operator itself is non-local), with q arbitrary, and the 
other functions equal to one. This Dirac operator has an eigenvalue spectrum on two parallel 
lines, the physical modes along the imaginary axis, and the doublers at infinite mass: mirroring 
the properties of the continuum operator. With the ability to use the function q to arbitrarily place 
the eigenvalues along the imaginary axis (except for the zero modes), this offers the tantalising 
possibility to construct a lattice Dirac operator whose eigenvalue spectrum precisely matches that 
of the continuum (provided that the lattice spacing is sufficiently fine that the lattice topological 
index matches the continuum topological charge). Since almost every observable in lattice QCD 
can be expressed in terms of the overlap operator and hence its eigenvalue spectrum, this 
operator would provide an intriguing possibility to reduce the artefacts of a continuum simulation: 
it could be a candidate perfect action. Of course, finding the crucial function q remained and 
remains elusive; and until it is found the hypothesis that such a perfect action can be constructed in 
this way remains questionable. 

Assuming that such a perfect action exists, a question arises concerning whether there is a con- 
nection between this purported 'perfect action' and the long-established perfect action derived from 
renormalization group considerations [Qj. This naturally leads to an examination of the place of 
overlap fermions within the language of the renormalization group. In this very preliminary study, I 
start an investigation of the place of overlap-style lattice Dirac operators within the renormalization 
group, in the hope of better understanding their theoretical foundation. These proceedings should 
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thus only be percieved as a few brief ruminations while I leave a more detailed construction to a 
future work JTot] - In section ||, I outline the required renormalization group tools; in section ||| I 
apply these to the case of overlap fermions on the lattice; and in section |] I use these results to give 



an alternative derivation of equation fll.lj). 



2. Block renormalization transformations 

I can define a block transformation, used to transform a spinor Yo an d Dirac operator Do to a 
spinor \\f\ and Dirac operator D\ in terms of three functions a, B, and B d : 

Z= d\if dyoe g o =N(a,p) d^d^i cfij7 d\j/Qe «r M x 

-(? 1 W-? (/)^ 1 (^^))«M(viW-B- 1 (yy)vo(y)) = f ^^j^ 1 ' 1 ' ^ 1 ' 1 '"?^, (2.1) 



where Af is a normalisation factor, and I will always assume integration (in the continuum) or 
summation (on the lattice) over repeated indices. I will also assume that B, B d and a are all 
analytic and invertible (in some sense; precisely what I mean by 'invertible' lies beyond the scope 
of this work). The derivation of the last equality only requires that B^aB^ 1 is positive-definite, 
that B d l aB~ Y +Do is invertible and that the 'inverse' blockings B and B d exist. 

Suppose that the original fermion action, \jf DoYo is invariant under an infinitesimal chiral 
rotation: 

Fo^vV* 75 ; Yb ( 2 - 2 ) 

Following the methods of Ginsparg and Wilson [p and Luscher [||], it can easily be shown that this 
implies a symmetry of the transformed action, where the new Dirac operator obeys the modified 
Ginsparg-Wilson relation 

D x B- x y 5 B + B d j 5 B^ l D x = D { {a- l B d y 5 B d l + B- l y 5 Ba- l )D u (2.3) 

and the corresponding 'chiral symmetry' transformations of the spinor fields are 

Y —>. e ie V>(.y5 B ~ 1 V5 B -V> B ~ 1 V> Ba ~ lD i)\i/ xjf ^T^g<£(Srf75Brf Vs-Oior'BrfrsB^'rs^ (2.4) 

while the topological charge is 

Q f = l -Tr{B- l y 5 B + B d y 5 B d l -B- l y 5 Ba- { D-Da- l B d y 5 B d l ) (2.5) 

If B and B d commute with 75, which both Ginsparg and Wilson's original work [jl]| and the con- 
struction of the standard forms of the fixed point action ^ assume then, unsurprisingly, equa- 



tion ( |2.3| ) reduces to the canonical Ginsparg-Wilson relations. Here I consider the different case 



where a is proportional to the unit matrix and has a value of 00 and B does not commute with 75. 



When a = °°, the integral over Yo in equation (2.1) can be performed explicitly, giving a 'Jaco 



bian' e~ J = det(Bj 1 B _1 ) (together with an unimportant constant factor) and a new Dirac operator 
D\ = B d DoB. Of course, on the lattice, for this chiral symmetry to reduce to the correct contin- 
uum expression, y^By=,B 1 and a should exist, be local, invertible, and reduce to 1 + 0(a) and 00 
respectively. 
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It is possible to construct a blocking transformation from the eigenvectors of the Dirac opera- 
tor. I introduce a large negative mass A, which will act as a momentum cut-off. For the continuum 
Dirac operator, the eigenvalue equations (for each pair of non-zero eigenvalues) can be written as 



(D o + A)0+(x,A) =(A + iX)(j> + (x,X) 
(D o + A)0_(x,A) =(A-iA)0_(x,A) 
0+(x,A) = y 5 0_(x,A) 



(2.6) 
(2.7) 
(2.8) 



From these eigenvalues and eigenfunctions, we can construct the eigenfunctions and eigenvalues 
of the Hermitian Dirac operator 75 (Do + A). The eigenfunctions are 



with eigenvalues ±jU where 

H=y/A 2 +X 2 



,i77(A,A) _ i?i+A 



VA 2 + A 2 ' 

We can use this formulation to construct a renormalization group blocking, with a = 00 and 



(2.9) 



(2.10) 



B 



B d = £ (j> (x , ,0)^ + fdXp{X) 

zero modes 



(j>-(x,X)(j>l(x',X)H 




+ 



(2.11) 



where p is the density of eigenvalues and H is an otherwise arbitrary analytic function satisfying 
3t(H[z} 2 ) > for 9t[z] > 0, which, given that 9t((l + e irt )/X) > ensures that the condition that 
the real part of the eigenvalues of B^ l B 1 must be positive is satisfied. For example, taking H = 1 
gives a new continuum Dirac operator 



Di = l + y 5 sign(y 5 (Do + A)), 



(2.12) 



where the eigenvalues of this continuum Dirac operator D\ lie on a semi-circle of radius 1 in the 
complex plane. The Jacobian for this blocking is 



f p(A)log 



1 



A 



2A 2 



VA 2 + A 2 A 2 _ 



(2.13) 



Assuming that we can regulate the Dirac operator so that we can neglect those eigenvalues which 
are not significantly smaller than A, we can expand J in powers of DqD/A, which gives 



/ = Tr 



tDI,D, 



29 



29 
8A4 



£F 2 v (x) + 0(A- 



(2.14) 



When I apply these methods to the lattice, assuming that Do is an ultra-local Dirac operator, 
Tr (D^D^) will be a constant (independent of the gauge field), and this Jacobian just gives a 
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j8— shift, as implied in equation (2.1). In the continuum, these trace terms can be re-expressed 
in terms of a bosonic field. 

Of course, D\ is only exponentially local, with a rate of decay controlled by A, and if we insist 
that in the continuum the Dirac operator must be 'ultra-local,' then this transformation would only 
be valid in the limit that A — > °°, where we recover the original operator Do- On the lattice, this 
situation is tolerable and controlable as long as A is inversely proportional to the lattice spacing. 
This exponential locality will always be seen if H is analytic and A negative [[Kj. This implies that 
the eigenvalues of the new Dirac operator will lie on a closed loop in the complex plane. In the 
general case, it is clear that if Do, Bj and B are local, then Di = BjDqB will also be local. 



3. The renormalization group and overlap fermions 



To move from the continuum to the lattice, one needs a renormalization group blocking which 
averages (in some respect) over the continuum fields. I introduce a blocking procedure P[(l,x), 
a function of the continuum gauge fields, so that the lattice fermion fields can be generated 
from the continuum fields Yc{ x ) by the relation YlQ) = / d 4 xP£(l,x)\j/c{x). It is also desirable to 
have an inverse blocking Yc(x) = E/^l 0M)Vk(O> so *b at the continuum theory can be extracted 
from the lattice theory by a continuous renormalization group transformation. These blockings 
have to be local, and include a means of removing the doublers in the corresponding Dirac operator; 
but otherwise we are very free in our choice of the blocking. Finding a blocking to move from the 
continuum to a particular lattice discretization is straight-forward. Finding a suitable inverse, to 
reach the continuum from the lattice, is not. For example, it is 'possible' to generate a Wilson-type 
fermion action from continuum fermions using the blocking 

¥c(y)=j:U d (\h-yp\-^) fi+Ir«</a-^)V" J ' AM(m VL(0=I^(y,0%(0 

x R V Z/ V a / x 



i-£y««(/a -y a ))e 



where n is the function 



n(x) 



1 *> (l-|e|)a/2 
\x\ < (1 — |e])a/2 
-1 x < (\-\e\)a/2 



i 

(3.1) 
(3.2) 



in the limit that e — > 0. Inserting ( p.l[ ) into the continuum action, \jr c (y')Dc(y' ,y)y/"c(y), gives 
an action ~\f/ L (l')D w (l' ,l)\f/ L (l). The only difference between this action and the standard Wilson 
fermion action is how the link U is constructed from the the gauge field : rather than a single 
parallel transporter, this new 'link' is a sum of parallel transporters. These blockings can be used 
as a basis for the transformations needed to generate overlap fermions. As an example, 

Bo,d =B w ,d', B =D w l (l + y 5 ^w^\gn{X w ))B w , (3.3) 

where (j>w and Aw are the eigenvectors and eigenvalues of 75 (Dw +A), generates the overlap action. 
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However, the Wilson-fermion 'blocking' given in equation (3. 1 ) is not valid, because although 



BjDcB gives an object similar to the Wilson operator, the matrix fi^S^ 1 is not positive definite; 
therefore, as mentioned earlier, the last equality of equation ( |1.1[ ) is invalid. However, for overlap 
fermions B^^Bq 1 is positive definite (provided that the lattice spacing is sufficiently fine); thus this 
is a valid renormalization group transformation. I shall leave the formal proof of this statement 
to the subsequent publication JToj] ; for now I will simply say that it is strongly suggested by the 
well-known fact that overlap fermions satisfy the Ginsparg-Wilson relation while Wilson fermions 
do not. 

There are alternative formalisms which can be used to construct the a and B blocking functions 
to obtain a particular lattice Dirac operator Dl, given some transformation P£ which averages 
over the continuum fields in a suitable manner to obtain the lattice fields. One possibility is the 
symmetric choice Bj = D L PlD c , B = D c PlD l > an d 05 = °°, which leads to a Ginsparg- 
Wilson relation 

= D L y 5 %y 5 + %D Ll (3.4) 

where 

7s = ^— 74- (3-5) 

d[d l 



It is easy to demonstrate that equation 3.4 is satisfied if we replace Di with the overlap operator. 
The Jacobian of the blocking is 



7 = Tr 



log(D c )-logp L )-log(P£F£) 



\a 4 F^ v + ..., (3.6) 



and again this transformation preserves the action, just giving a shift in the coupling constant. This 



choice of B is not unique Qll|]; but the work here suggests that the different ways of expressing 

the chiral symmetry on the lattice are related to different lattice regularisations (since each of these 

possible blockings will give a different shift to the gauge coupling). For example, if we choose 
— C — i r 

Bd = P L , B = D c P^Di, we recover the canonical Ginsparg-Wilson relation. It is to be stressed 
that unlike Ginsparg and Wilson's original chiral fermions, and perfect action fermions which use 
a different blocking with a finite a, the lattice chiral fermions discussed here are generated with 
a = oo and a blocking field B which does not commute with 75. 

4. New Solutions to the Ginsparg-Wilson operator 

Once we have found one chiral lattice Dirac operator, it is easy to generate more by construct- 
ing blockings from the eigenvalues of the operator. For example, by witing the overlap eigenvalues 
as 1 + e' 9 with eigenvectors 0, we can construct a blocking 

Ydy) =^F(G)&(d,y)&He,x)x l / L (x) 

x e 

WM =I? , LWlf(e)®V^)0(ej) 



1 +e w s(e 



' ?(cos0)%(cose)(l+r(cos0y e )] (4-1) 
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where s,t,h, q and r are analytic functions (chosen to ensure the correct continuum limit and that 
the fermion doublers have infinite mass in the continuum, which can be achieved by ensuring that F 
is analytic, positive, and that F = 1 at 6 = 0, 7i,2n), which leads to the generalised overlap operator 



given in equation That this operator is exponentially local is clear from the analyticity of the 
momentum representation of the renormalization group transformation, the exponential locality of 
the overlap operator and that the overlap operator commutes with its Hermitian conjugate. 

5. Conclusions 

I have tentatively suggested that there is, in an infinite volume, a renormalization group trans- 
formation linking the continuum fermion action with a Ginsparg-Wilson lattice fermion action, and 
I have shown how such a blocking can be constructed for massless overlap fermions. However, a 
different style of block transformation is used compared to that in fixed point fermions or original 
construction of the Ginsparg-Wilson relation. The Jacobian from the block transformation corre- 
sponds to a shift in the gauge coupling, together with some irrelevant terms; thus the lattice QCD 
Lagrangian is preserved by this transformation. 

Clearly there is one crucial step missing in the argument as outlined here, namely a proof that 
for the overlap operator the matrix, defined here as B^ y B l , is positive definite. I shall attempt 
to provide that proof in a future publication However, if this conclusion is correct, then it 
could have implications in the analysis of discretization errors, and in the theoretical description of 
a perfect action derived from overlap fermions. 
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